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Abstract 

The many-body wave scattering problem is studied in the case where the bodies 
are small, ka <C 1, where a is the characteristic size of a body. The limiting case when 
a — > and the total number of the small bodies is M = 0(a 2 ~ K ) is studied. 
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1 Introduction 

The many-body scattering problem in the case of small scatterers embedded in an inho- 
mogeneous medium has been solved in [1] and [2] under the following assumptions: 

ka«l, d = 0{a^), U = f ^, (1) 

where a is the characteristic size of the small bodies, n = 2tt/X = uj/cq is the wave number 
and Co is the wave speed in free space, k 6 (0, 1) is a parameter one can choose as one 
wishes, d is the distance between neighboring particles, h(x) is a piecewise-continuous 
function in a bounded domain D C M 3 with a smooth boundary S, = hi < 0, h = 
hi + i/12, x m G D m is an arbitrary point, D m is a small body, S m is its surface, N is the 
unit normal to S m , 1 < m < M, M is the total number of the embedded small bodies in 
D, the unit normal N points out of D m , Cm is the boundary impedance in the boundary 
condition: 

du 

— = CmU on S m , 1 < m < M; u = u M , (2) 
and the distribution of small bodies in D is defined as 

AA(A) := S DmCA l = ^ / iV(x)dx[l + o(l)], a -> 0, (3) 
or K J a 
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where A/"(A) is the number of small bodies in an arbitrary subdomain Acfl, -^( x ) > 
is a piecewise-continuous function, and for simplicity it is assumed that D m = B(x m ,a) 
is a ball centred at the point x m , of radius a. The scattering problem, solved in p]-[2], 
consisted of finding the solution to the equation 



M 



[V 2 + k 2 n 2 {x)]u m = in M 3 \ (J D m , (4) 

m=l 

satisfying boundary conditions ([2j) and the radiation condition 

Um = u + v m , ^t^- ~ ikv M = o(-), r := |x| -> oo. (5) 
or r 

Here, no is the solution to problem (HJ)-© in the absence of the embedded particles, 
i.e., the solution for the problem with M = 0: 

[V 2 + k 2 nl(x)]u = in R 3 , (6) 

where n^{x) is the refraction coefficient in the absence of embedded particles, tlq(x) = 1 
in the region D' := M 3 \ D, and 

n = e ifcQ - x + v , ^-ikv = o(^], r -» oo, (7) 

where a G 5* 2 is the direction of propagation of the incident plane wave, and S 2 is the unit 
sphere in R 3 . 

It was proved in PQ, [2], that, as a — > 0, the limiting field n does exist and solves the 
equation 

[V 2 + A;V(x)]u = in M 3 , (8) 

where 

n 2 (x) = n§(x) - 4vr/c~ 2 /i(x)iV(x). (9) 

Therefore, in the limit a — > 0, under the constraints - ([3]), the limiting medium, 
obtained by the embedding of many small particles, has the refraction coefficient n 2 (x), 
given by ([9]). Since the functions /i(x) and iV(x) are at our disposal, subject to the 
restrictions iV(x) > 0, 9/i(x) < 0, it is possible to create any desired refraction coefficient 
n 2 (x), 3 : n 2 (x) > 0, by choosing /i(x) and iV(x) suitably. 

It is assumed that the term "piecewise-continuous" function / in this paper means 
that the set M. of discontinuities of / is of Lebesgue's measure zero and, if S is a subset 
of this set such that / is unbounded on S, f\s = oo, then / grows not too fast as x tends 
to S: 

l/(x)l * [dist^r °- l/<3 ' c=const ^ ' ( 10 ) 

so that the integral f D /(x)<ix exists as an improper integral. 

This paper is closely related to PQ, and its goal is to develop a theory, similar to the 
one in [I], for a different governing equation: 

LqUq := V • (c 2 (x)Vu ) + uj 2 u = in R 3 , (11) 
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where the wave speed c(x) = Co = const in D' := R S \D, the complement of D in i? 3 , and 
c(x) is a smooth and strictly positive function in D. The speed c(x), in general, has 5" 
as its discontinuity surface. In this case equation (jlip is understood in the distributional 
sense as an integral identity: 

/ (-c 2 (x)V0Vu o + w 2 <M))dx = G C^^ 3 ). (12) 

Alternatively, one may understand equation (jlip as the following transmission problem: 

LqUq = in D, = uo in D, (13) 
LoUq~ = in D' , Uq = uq in D' (14) 

U+ = Uq, c 2( x )^. = c 2( x )^o_ on S . (15) 

The transmission conditions (115j) together with equations fjl3[> and (I14p are equivalent 
to problem (fT2|) . Existence and uniqueness of the solution to (fl~3]) - (fT5j) was proved in 

The scattering problem we are interested in can be stated as follows: 

M 

L u = in R 3 \ [J D m ; u = u M , (16) 

m=l 

— = ( m u on S m , 1 < m < M, (17) 

u = uo + v, v r — ikv = o(-), r — > oo. (18) 

r 

In section [2] problem (|16p - (|18|) is investigated and the limiting behavior of u as a — > 
is found. 

We conclude this Introduction by a brief derivation of the governing equation (fTTj) . 
The starting point is the Euler equation: 

v + (v,V)v = -^, (19) 
P 

where v is the velocity vector of the sound wave, p = p(p) is the static pressure, p is the 
density, and 

Vp = c 2 (x)Vp, (20) 

where c(x) is the sound speed. 

Let the material in D be initially at rest, v = v(x, t) be a small perturbation of the 
equilibrium zero velocity, the density be of the form p = po + VK x ^)i where po is the 
equilibrium density of the material, which is assumed to be constant, and ip and v are 
small quantities of the same order of smallness. 

The continuity equation is 

^ = -V • (pov), (21) 
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where the term V • (iftv) of the higher order of smallness is neglected. Differentiating (|2ip 
with respect to time yields 

4> = -V • (pov). (22) 

Under the same assumptions about p = po + ^(x, t) and v, the term (v, V)v in (|19p 
is of the higher order of smallness and is therefore neglected. Multiplying (|19p by p and 
neglecting the term ipv of higher order of smallness yields the acoustic momentum equation 

p v = -Vp. (23) 

Substituting ([2D ]) in ([25]1 gives 

p v = -c 2 (x)V^, (24) 

where the relation Vp = was used. This relation is exact for a constant po- 
Substituting (J2JJ) in ([22J yields 

^- V- (c 2 (x)VV>) = 0. (25) 

If -0 = e~ tult u, then ([25]) reduces to equation (ITT]) . 



2 The scattering problem 

In this Section, problem (fTB"]) - ([18]) is studied. Assumptions ([JJ and ([3]) are still valid. 
Let G be the Green's function for the operator Lq: 

L G(x,y) = -<5(x-y) in R 3 , (26) 

G satisfies the radiation condition 

dG 1 

—. — r — ikG = of: — r ) as Ixl — > oo. (27) 
a|x| |x| 

The following result from [3] will be used. 

Theorem 1. In a neighborhood of a point of smoothness o/c(x) one has 

G(x,y)= 1 r (l + o(l)), |x-yH0. (28) 

47r|x — y|c(x) 

In a neighborhood of the point x £ S 1 , where S us a smooth discontinuity surface of 
c(x), one /ias 

f rri J tvJ + +o(l)], ye D, 
Gxj) = ^ 4 ™+W! x ^ * ' (29) 

li^y^xy 1 -^ 1 + o(i)], ye i?'. 1 ; 



where 



6:= c + (x) c-(x) rxy;=|x _ y|) i?= Vp 2 + (|x 3 | + |y 3 |) 2 , (30) 



V^i-yi) 2 + (x 2 -y 2 ) 2 . (31) 
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The origin of the local coordinate system lies on S, the plane X3 = is tangent to S, 
c+(x) and c_(x) are the limiting values o/c(x) when x —* S from inside (and outside) of 
D. 

In [3], the operator Lq corresponds to the case u = 0. However, in p] it is proved that 
adding to Lq a term q(x)G(x,y) with a bounded function q does not change the main 
term of the asymptotic of G as x — » y. 

The solution to problem (I16j) - (1151) is sought in the form 



M 

U = Uq 



M r 

+ Y. G(x,t)a m (t)dt. (32) 

m=l Sm 



For any a m G L 2 (S m ), the function u, defined in ([32]) . solves equation ([16]) and satisfies 
the radiation condition (|15p , since G does. Therefore, (|32p will be solution to problem 
(fl"6|) - (fT5|) if a m are such that the boundary conditions (fTT|) are satisfied. Uniqueness of 
the solution to problem (116j) - ()15f) follows from essentially the same arguments as in PQ, 
see the proof of Theorem 1 in [1] . 

The boundary conditions (|15j) imply 



where 
and 



UeN + Am(Tm - °" mCm = C m 1i e + Cm / G(s, t)(T m (t)dt, (33) 
• Cm ^(^m)/& j 

ue(x) := := «o(x) + V / G(x, t)<7 ro / (t)tft. (34) 

— ./ 9' 



The field iti" 1 ^ is called the effective (self-consistent) field. It is the field acting on the m-th 
particle from all other particles and from the incident field uq. 

The operator A m is the operator of the normal derivative of the single-layer potential 



Ta m := [ G(x,t)a m (t)dt 

J Sm 



at the boundary S, and 
dTa m Aa m 

) „ r dG(s,t) 



Aa m = j } ] Pm{t)dt, seS ra , (35) 



In equation (|35p , is the limiting value of the normal derivative on S m from outside 
of D m . Equation (135f) is well known from the potential theory in the case c(x) = 1 and 

G ( x >y) = r xy : = l x -yl- 

If c(x) 1, then by Theorem 1 one may consider T and A as the operators, cor- 
responding to c(x) = 1, divided by c(x m ), because c(s) is assumed smooth in D, and, 
therefore, it varies negligibly on the small distances of the order a. 
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The basic idea of solving the many-body scattering problem (|16p - f)18|) is similar to 
the idea originally used in pQ. 

The approach is to reduce the solution of the many-body scattering problem by small 
bodies to finding some numbers, rather than the unknown functions a m , 1 < m < M. If 
M is very large, it is practically impossible to use the usual system of boundary integral 
equations for finding the unknown a m . 

let us rewrite equation ([32]) as follows: 



where 



G(x, x m )Q m + / [C(x, t) - G(x, Xm )]a m (t)dt, (36) 

m=l m=l Sm 

Qm ■= / a m (t)dt, (37) 



and prove, as it is done in [T], that 



G(x, x m )G/ r 



[G(x,t) -G(x,x m )]a m (i)dt 



a — > 0; |x — x m | > a. (38) 



Thus, the solution u can be written as 

M 



u = Uq 



(x) + G(x, x m )Q m , |x — x m | > a, 



(39) 



m=l 



with the error that tends to zero as a — » 0. 

Consequently, the scattering problem is solved if the numbers Q m , 1 < m < M, 
are found. This simplifies the solution of the many-body scattering problem drastically, 
because equation ([32]) requires the knowledge of the functions cr m (i), 1 < in < M, rather 
than the numbers Q m , in order to find the solution u of the scattering problem. 

The next step is to derive the main term of the asymptotics of Q m as a — > 0. 

To do this, we integrate (|33p over S m and neglect the terms of the higher order of 
smallness as a — > 0. 

One has: 



/ u eN ds= / V 2 u e dx = (V 2 u e )(x m )|D m |, \D r , 

JS m JD 

ill 



47ra 3 



(40) 



where the Gauss divergence theorem was applied and a mean value formula for the integral 
over D m was used. 
Furthemore, 

Qm 

Cm 



1 



/ Aa m ds = -—j a m ds 



(41) 



where (cf [T] 



/ ds [ d 1 a U)dt = — - f a (t)dt. 

JS m JSm dN * 4vrr st c m Js m 



Aads :=-?-/ ds . . ( _ 



(42) 



6 



Thus, integrating ([33]) over S m yields 

V 2 u e (x m )|L> m | - c m x Q m = ( m u e (x m )\S m \ + — f dta m (t)[ ds—^—, (43) 

Cm Js m JS m 4TTr st 

where \S m \ = Ana 2 is the surface area of the sphere S m and formula ([28]) was used, namely, 
we have replaced G(s, t) by 47r 1 r - using the smallness of D m , and we have replaced c(s) 
by c(x m ) = c m because \x m — s\ < a, and a is small. 
Using the identity 

ds 

a if |s — x m | = a and \t — x m | = a, (44) 



one gets from ([43 p the following relation: 

Qm{c m l + c^Cmo) = -47rC m u e (x m )a 2 + 0(a 3 ). (45) 

If a — > and k £ (0, 1), then 

( m a = h(x m )a 1 ~ K = o(l), a ->• 0, 

the term (9(a 3 ) in (|45jl can be neglected, and one gets the main term of the asymptotics 
of Q m as a — > 0, namely: 

Q m = -4vr/i(x m )n e (x m )c(x m )a 2 - K [l + o(l)], a -> 0. (46) 

Therefore ([34]) , ([39]) and ([46]) yield 



u e (x) = u (x) - 47r G(x,x m )/i(x m /)u e (x m /)c(x m /)a 2 K [l + o(l)]. (47) 

m'^m 

Taking x = x m and neglecting o(l) term in (|47p . one gets a linear algebraic system for 
the unknown quantities u m := u e (x m ), 1 < m < M , 

u m = u 0m ~ 47T ^2 G(x m ,x m /)/i(x m /)c(x m /)n m /a 2_K . (48) 

Let us now derive and use a generalization of the result proved originally in pQ. This 
generalization is formulated as Theorem 2 below. 
Consider the sum 



M 

/=lima 2 - K V/(x m ), (49) 

a^O 1 — » 

m=l 

where the points x m are distributed in D according to ([3]) . 

Assume that /(x) is piecewise-continuous and (jlOp holds. If / is unbounded, that is, 
the set S is not empty, then the sum ([4^]) is understood as follows: 



M 

I := lim lim a 2 ~ K V /(x m ). (50) 

<5— »0 ' — • 

m=l,dist(x m ,cS)>5 
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Theorem 2. Under the above assumptions, there exists the limit (|49p and 



M . 
lim a 2 - V /(x m ) = / /(x)iV(x)dx. 



(51) 



Proof of Theorem 2 is given at the end of this paper. 

Applying Theorem 2 to the sum (|47p one obtains the following result: 

Theorem 3. There exists the limit: 

lim u e (x) := u(x), 

and the limiting function solves the equation: 

u(x) = u (x) - 4tt / G(x,y)h(y)c(y)N(y)u(y)dy. (52) 

Applying operator Lq, defined in (fTTj) , to ([52]) and using the relations 

L G = -S(x-y), L u = 0, (53) 

one obtains the following new equation for the limiting effective field u: 

L u = 4vr/i(x)c(x)A^(x)n. (54) 

This equation can be written as: 

Lu := V • (c 2 (x)Vu) + u?u - 4vr/i(x)c(x)A^(x)ti = 0. (55) 

Therefore, embedding many small particles into D and assuming (H|) - © , one obtains 
in the limit a — > a medium with essentially different properties described by the new 
equation 



Let us now prove Theorem 2. 
Proof of Theorem 2. 

Let S be the subset of the set of discontinuities of / on which / is unbounded, let the 
assumption (|10p hold, and let 

D s :={x:xe D, dist{x, S) > 5}. (56) 

Consider a partition of Dg into a union of small cubes A p , centered at the points y p , 
with the side b = a 1 / 3 . One has 

M 

/(x m )=]T/(y p )[l + (l)]a 2 - K £ 1 
m=l,dist(x m ,M)>8 P x m €A p 



J2f(y P )N(y p )\A p \[l + o(l)} (57) 
p 

> I f(y)N(y)dy as a -> 0. 
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Here in the second sum we replaced /(x m ) by f(y p ) for all points x m 6 A p . This 
is done with the error o(l) as a — * 0, because / is continuous in D$. In the third sum 
we have used formula ([3]) for A = A p . The last conclusion, namely, the existence of 
the limit as a —* 0, follows from the known result: the Riemannian sum of a piecewise- 
continuous bounded in D$ function /(x)iV(x) converges to the integral J D f(x)N(x)dx 
if maXpdiamAp — ► 0. In our case 

diamAp = V3a 1/3 as a -> 0, (58) 

so formula (|57p follows. 

From the assumption ()10|) with v < 3 one concludes that 

lim / f(x)N(x)dx = [ f(x)N(x)dx. (59) 

The integral on the right in (|59p exists as an improper integral if v is less than the 
dimension of the space, i.e., v < 3. Therefore, formula (|5ip is established. 

Theorem 2 is proved. □ 
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